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We generalize the Landau-Khalatnikov hydrodynamic theory for superfluid helium to two- 
component (binary) Bose mixtures at arbitrary temperatures. In particular, we include the spin-drag 
terms that correspond to viscous coupling between the clouds. Therefore, our theory not only de¬ 
scribes the usual collective modes of the individual components, e.g., first and second sound, but 
also results in new collective modes, where both constituents participate. We study these modes 
in detail and present their dispersions using thermodynamic quantities obtained within the Popov 
approximation. 
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I. INTRODUCTION 

The realization of a Bose-Einstein condensate (BEC) 
in ultracold alkali-metal vapors m has ignited a rapid 
progress in the understanding of degenerate gases at low 
temperatures HH]. A large part of this understanding 
has been gained through the study of collective modes 
[6H8]. In particular, it has been shown that for a weakly 
interacting gas of bosons close to absolute zero tempera¬ 
ture, the collective excitations are Bogoliubov quasiparti¬ 
cles [sna that are responsible for fascinating properties 
of the system, including superfluidity and quantum de¬ 
pletion of the condensate. 

Work building upon the single-component ultracold 
gases has provided access to even richer systems. In par¬ 
ticular, considering mixtures of several species of parti¬ 
cles with the same (Bose-Bose or Fermi-Eermi mixtures 
[inmi]) or different statistics (Bose-Fermi mixtures [15]) 
has become possible. These systems are known as binary 
mixtures or two-component gases. Arguably the simplest 
of them is a mixture of two different hyperfine states of 
the same bosonic atom. However, even this simple sys¬ 
tem poses important questions concerning the nature of 
its ground state and the excitations. Therefore, much 
work has been carried out on the static and dynamic 
properties of the two-component Bose gas [HHZS], both 
in the uniform case and for the trapped case. Most of 
the effort has been concentrated on the zero-temperature 
physics, with only a few studies (SOI ED on the properties 
of binary Bose mixtures at nonzero temperature. 

Having more than a single component in the gas also 
allows one to make a connection to the physics of spins, 
by introducing a pseudospin to distinguish the two com¬ 
ponents. In particular, one can consider ferromagnetic 
and antiferromagnetic states as well as spin dy¬ 

namics [35l |36] , and topological spin textures [32l |33l l37l - 
l44] . One kinetic effect concerning the spin dynamics is 
the so-called spin drag [151135] . This recently observed 
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SZlllHl effect corresponds to the relaxation of the differ¬ 
ence of the velocities between the two components. Un¬ 
derstanding the interplay between the BEC, the thermal 
particles, and the spin degrees of freedom in this rel¬ 
atively simple and well-controlled two-component Bose 
gas might also offer some insights for interacting spinful 
degenerate systems of a rather different nature, such as 
the condensate of magnons |49l |50] . 

In this paper we build upon our previous results for 
the ferromagnetic Bose gas [22], but now consider a dif¬ 
ferent situation where two condensates are present in the 
miscible regime. We tackle the problem of the collec¬ 
tive modes of the two-component mixture both in the 
uniform gas and in a trap in an effort toward making a 
connection with experiments. The structure of the paper 
is as follows. In Sec. H we describe the Popov theory 
of the binary Bose mixture, and present relevant ther¬ 
modynamic functions in that approximation, including 
the equation of state. We study the effects of spin drag 
in Sec. HI. In Sec. IV we develop a linear hydrodynamic 
model that makes use of our previous thermodynamic re¬ 
sults and describes a two-component system accounting 
for spin drag. We present the results for the uniform and 
trapped gas in Sec. V. Finally, we conclude in Sec. VI. 


II. MICROSCOPIC THEORY 

In this section we briefly describe the microscopic 
Popov theory of the two-component Bose mixture. The 
Popov theory is an extension of the Bogoliubov theory 
to relatively high temperatures, which includes an im¬ 
proved treatment of the excitations. Specihcally, the Bo¬ 
goliubov excitations are allowed to interact, and their 
interactions are treated in the Hartree-Fock approxima¬ 
tion. Multi-component gases of bosons have been treated 
in the Bogoliubov framework before [51] . In particular, 
two-component mixtures have been considered in Refs. 
[HUM ED, and some results from the Popov theory 
have been presented in Ref. m- The novelty of our 
results is twofold: We present the Popov analysis in the 
functional-integral formalism and calculate the thermo- 
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dynamic properties of the balanced binary Bose gas. Our 
discussion on the Bogoliubov transformation follows the 
usual grand-canonical treatment of the problem. The 
single-component situation has been treated in this way 
in, for instance, Refs. [551 [M] . 


A. General binary mixture 

In general, a grand-canonical partition function for two 
bosonic fields (0^ and that includes all the possible 
s-wave interactions is 


where the action is 

rhj} p 

4)^, 4)1, (/);] = J dr J dx 


( 1 ) 

( 2 ) 
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and all the fields are considered at the position x and the 
imaginary time r. Moreover, /? = l/k^T is the inverse 
thermal energy, Wa are the masses of the particles, and 
ga-a' are the two-body T matrices describing the s-wave 
interactions. 

From now on we focus on the symmetric case, where 
the masses are equal m = = m^, and the intraspecies 

interactions are the same and described by a single scat¬ 
tering length a. Thus, 

Anfisa 

9 = 9n=9u = -• ( 3 ) 

m 

The interspecies interactions are described by another 
scattering length implying that 


9U - 


m 


(4) 


All the interactions are assumed to be repulsive, i.e., g > 
0 and g^i > 0. Furthermore, as opposed to our earlier 
work in Ref. 1121, we here focus on the case with two 


separate condensates. In order for this to be possible, 
the condition 


9 > 9\i 


(5) 


has to be satisfied, as otherwise the two components 
demix [55] . 

We are now in a position to perform a fluctuation ex¬ 
pansion for each species by putting 

4),y{X,T)=4)0,y{x)+4)'^{X,T), ( 6 ) 

where the fluctuations 4 ’'^{x,t) are, on average, zero. 
Moreover, the fluctuations are orthogonal to the conden¬ 
sate (0a(®, t)) = 0oa(®) of the same species which means 
that 

J dx(^4)*Q^{x)4)'^{x,T) + 4)Oa{x)4)'*{x,T)^ =0. (7) 

Since in what follows the relative phases of the conden¬ 
sates do not play a significant role, we choose both the 
condensate fields to be real, 

0Oa(^) — ^Oa; (8) 

where noa is the atomic (number) density of the con¬ 
densed a particles. Moreover, since we are considering 
the uniform case here, the condensate density has no spa¬ 
tial dependence. Expanding the fields in the action in 
this manner, we have for the action 

‘^ = '^0+ ^ 5ia + 52 + 53 + ^4, (9) 

a=t.t 

where the zeroth-order (Landau-free-energy) contribu¬ 
tion is 


So = hjdV ( - -k 

- + 5’^ot/2 + , (10) 

the term linear in fluctuations reads 

nhl5 j. 

5la[0a>0a]=y drjdx 


0a {- 9 a + gnoo 


+gn^os)V^ + c.c. 


( 11 ) 


and the quadratic term is 


52[0(.%0^,0(*,01] 




dx 


E 

a=t,t 
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^a + 2gnoa + gu'n-oa^ 0a 


g^Qq 
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+ guV^otnoi J dr J dx {4)\4)[ + 4)'^4)[ + 4)\4)'* + 4)'^4 )'^), 


( 12 ) 


where all the fluctuation fields are evaluated at {x,t), Furthermore, S 3 and 54 terms describe the interactions 

and we have denoted the species opposite to a by tt. between the fluctuations, 
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= / dr fdxf g^nocr {(j)'*<!)'* 

do J \ 


> ( 13 ) 


and 



Note that the terms S 3 and 5*4 are neglected in the Bo- 
goliubov theory. 

We now perform the Hartree-Fock theory for the ex¬ 
citations, which involves the inclusion of the mean field, 


(15) 

where n'^ is the density of the excitations of the a species 
such that the total density of a species is 

Ucr = nocr + n'^- (16) 

Note that we neglect the coherence between the two 
species and take {(j)'*= 0. Therefore, the appropriate 
mean-field substitutions are 

^ - 2n;^ (17) 

<(’t ~ (13) 

where the subtractions account for double counting in 
the quartic term of the action, whereas in the cubic terms 
no double-counting problems appear, as can be seen by 
applying Wick’s theorem. 

By requiring all the linear terms in </)(, and (j)'* of the 
action to vanish, we obtain a set of two Gross-Pitaevskii 
equations for the uniform condensates that read 

{-dt + gno^ + g^noi + 2gn\ + = 0’ (19) 

{-di + gnoi + 5Tl«ot + ‘dgn\ -f 3n^)v^ = 0, (20) 

and from which the chemical potentials are obtained as 

+ "d-gn'^ + (21) 

gi = gnoi + 5 n”-ot + 25’^! + (22) 

In order to diagonalize the quadratic part of the action, 

we perform a Fourier transformation, and then introduce 
Nambu space [54] . Since we want to rewrite the quadratic 
part of the action in the form 

C- '(’1] = -^ E (23) 

fc^0,n 

where hk is the momentum, n labels the Matsubara fre¬ 
quencies LOn = 27rn/fi,/3, 

^kn = {(l>'*km ^'*km ^l-kn) i (24) 

is a vector in the appropriate Nambu space in this case, 
and G~^ is the inverse Green’s function of the system. 


r 


Note that we have to take care to preserve the correct 
time ordering. The latter results in an extra term in the 
action, 

hp 

^TO = - ^ E [^k- g<j+ gi2noa- -b 2n'^) 


(7,k^0,n 


+5n(”'0rT +0] =-^ E i^k+gnoa), (25) 

(T,k^0,n 

where Sk = /2m is the kinetic energy. Therefore, 

plugging the expressions for the chemical potentials 
into the action, we have 


S = - hpV (^gun^ni + g{r 


'I ■ 




(26) 


Tin I +n/ 


g 


‘0^ 


ot 


-V E 


' gnoa) 


-X E 


(T.fe^O 

k^0,n 

where V is the volume of the system, while 






HT, 


■n 




(27) 


where and are two-by-two submatrices (from 

now on two-by-two matrices are denoted by capital let¬ 
ters, while four-by-four matrices are denoted by bold cap¬ 
ital letters). The submatrices on the diagonal are exactly 
the same as the Bogoliubov (single-component) inverse 
Green’s functions, that is. 


-f^GYku 


-itujn + Sk+ gnoa gnoa 

griQ^ ihwn + ek+ gn^^ 


(28) 

and the off-diagonal matrix is the self-energy due to the 
interspecies coupling given by 

'l 1^ 


= gtiVdPynol 


1 1 


(29) 


Anticipating the Bogoliubov transformation we define 


/ = 


where 


az 0 

0 cr. 


1 0 
0 -1 


(30) 


(31) 


Moreover, we define a matrix which is obtained by 
setting ojn to zero in the inverse Green’s function, i.e., 


Ffe = -hG 


kO ■ 


(32) 
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We now proceed to diagonalize the matrix while keep¬ 
ing the bosonic character of the excitations, i.e., requir¬ 
ing that their operators obey bosonic commutation rela¬ 
tions. In our formalism the preservation of the bosonic 
commutation relations is represented by the fact that the 
Matsubara-frequency terms are unaffected by the trans¬ 
formation. Hence, we look for real eigenvectors w of 

I Tk, 




satisfying the property 

w±^k ■ I ■ w±^k = 1. (34) 

The above-mentioned eigenvalues define the dispersions 
of the Bogoliubov quasiparticles 

+ 5(«ot + ^0^)) (35) 

± £k\j9'^{no\ - 

The eigenvectors w± k have the so-called Bogoliubov co¬ 
herence factors as their entries, or explicitly 

w±,k = (M±,fc, (36) 

For the case at hand these coherence factors are rather 

involved functions of the momentum. Since they have 
been presented explicitly in Ref. [24) . we do not write 
them out here. 

The relation between the two original fluctuation helds 
and the new Bogoliubov quasiparticle helds is given by 

^kn = Wk-^kn, (37) 


where we have dehned the Bogoliubov helds 

^kn = ('0-|-,fen> V'-l-, —fcri) V'—,fcn) V'—, —fcn) ) 

and the transformation matrix 


lVk = 


w 


+ .k 


w 


+ ,k 




which consists of the submatrices 


S,k 1 _ yC 


-^S.k 


<,k K,k 


(38) 

(39) 


(40) 


It is straightforward to check that this transformation 
leaves the Matsubara-frequency terms in the action un¬ 
affected, and therefore the Bogoliubov excitations are 
bosons. The action becomes 

S = - npV + g{nl + n^) - ^ 


+ ^ {E+:k + — 

k^O 


2ek - g(not + no;)) 


+ E E 

s—± fc^0,n 


where (7?+.fc -I- E_ f,) in the second term is again due to 
the time ordering but this time of the '4’± helds. 

In the preceding equation, the hrst term describes the 
Hartree-Fock contribution to the action. The second 
term (after properly accounting for the fact that the con¬ 
tact potential does not fall off at high momenta) can be 


shown to describe the so-called Lee-Huang-Yang correc¬ 
tion [MUST], which is small for a weakly interacting gas, 
and can therefore be safely neglected. The last term de¬ 
scribes the Bogoliubov excitations. Evaluating this path 
integral amounts to a Gaussian integration and can be 
performed exactly. Finally, we perform the remaining 
bosonic Matsubara sum 


lim In [/3(-zIiw„-I-Es_fc)] 

77—>-0+ 

n 

= ln(I-e-'^^=-'=), (42) 

and arrive at the following expression for the partition 
function 


Z = exp 


/3(g;;n;n; -f g(n: 


2', 

■«t) - g — 


'■OtN 


1 

V 


2 

-PE, 


ln(I — e 

s—±,fc^0 


(43) 


The partition function is related to the pressure by 
p{no^, noi, T) = -Vl/V = ^ ln(Z) 


= -F g(nj -h n^) - g 




*0; 


~ E ln(l-e ^^“■'=), (44) 

P s=±,k^0 

where fl is the grand potential. We have obtained the 
(average) particle densities 


ria = Tiao + ^ 


|2 


s—±,k^0 


gPE, _ X 


(45) 


from the appropriate Green’s functions, while the entropy 
per volume [58] is 


5 

V 


dfi 

df 


ks 

Y 

ks 

V 


y^ ln(l — e 

s—ib,A;^0 


^0Es,k _ \ ' 

s=ib,fe^0 


(46) 


For future convenience, we define the entropy per particle 
as 


S 1 
Vn^ + ni 


(47) 


and also the total chemical potential ptot, as well as the 
difference of the chemical potentials 

gtot=gt + g;, (48) 

Ag = - g;. (49) 

Similarly, we define the total particle density and the 
difference of the particle densities: 

ntot=n^+ni, (50) 

An = n; — n;. (51) 

Finally, it is also beneficial to define two additional T 
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matrices, 


g±= g±g^l= g{l±-f), (52) 

where we have also introduced a dimensionless number 
7, which shows the relative strength between the inter¬ 
species and the intraspecies repulsion. Note that in the 
miscible case that we discuss here, 7 < 1 . 


B. Balanced mixture 


Of particular interest is the balanced case, where the 
number of particles of the two species are equal: = 

= n. We consider it in this subsection. An obvious 
consequence of this limit is 

«ot = = ^^o, (53) 

and 


/^ = = Mi = 9 +n + gn. (54) 


Moreover, the dispersions of the quasiparticles become 
= £fc(efc 2nog±), (55) 

and the Bogoliubov transformation matrix simplifies con¬ 
siderably to 


Wk = 

which can be inverted to 

51 


W+,k 

,^+,k 

W..k] 

w-.kj ’ 

(56) 

( W- 

w:^k) ’ 

(57) 


from which we can conclude that the ijj+,kn held has equal 
contributions from the (j)^kn and (j^ikn fields, whereas the 
latter fields enter ip-,kn with a relative minus sign but 
with equal absolute weights. This implies that ijj+.kn 
describes densitylike excitations, whereas V'-.fe describes 
spinlike excitations. The submatrices Ws^k in the above 
are 





-Vs,k\ 

‘^s.k j 


(58) 


where in this case the coherence factors are simple enough 
to be written out explicitly as 


- 4 
1 
4 




V E±,k 
'ep+g±r 

^ E±,k 


(59) 

(60) 


Note that the coherence factors are very similar to the 
single-species case. However, the prefactor (1/4) here 
is different from the single-species case ( 1 / 2 ), since the 
transformation now involves four fields instead of two. 

We now proceed to discuss the thermodynamic func¬ 
tions of the balanced binary mixture. Throughout the 
discussion, we consider three different dimensionless in¬ 
teraction parameters = 0.01, 0.05, and 0.1. They 
have been chosen to correspond to the experimentally 


relevant weakly interacting (far away from Feshbach res¬ 
onances) ultracold gas situations. In particular, we 
consider the sodium atom which has several scattering 
lengths between the accessible hyperfine levels close to 
50 Bohr radii m- Moreover, we are interested in the 
hydrodynamic regime, where the density in the center of 
the trap might become as high as 10 ^^ m’ which 

corresponds to ~ 0.03. Comparing thermodynamic 
quantities calculated within the Bogoliubov theory with 
the renormalization group results for the single-species 
case (see Ref. M) shows that the two agree very well 
all the way up to the condensation temperature Tc for 
= 0.01. Agreement for = 0.05 and par¬ 

ticularly for 0.1 is less good. However, the qualitative 
features of the thermodynamic functions are preserved. 
We expect the situation to be similar for the two-species 
case. 

We present the equation of state (pressure) in Fig. 
We scale the pressure plots with the zero-temperature 
pressure g+n^/2 obtained from a Gross-Pitaevskii cal¬ 
culation [55) . For increasing interaction strength, the 
scaled pressure decreases, since the thermal energy be¬ 
comes comparable with the interaction energy at higher 
temperature. For 7 = 0 , the pressure is equal to the 
pressure of a single species of gas with twice the density 
(see, e.g.. Ref. [SI] for comparison). 

Note that for a fixed total number of particles, the pres¬ 
sure is not a monotonically increasing function of temper¬ 
ature in the Bogoliubov theory, as opposed to the Popov 
theory. This spurious effect appears due to the competi¬ 
tion between the depletion of the condensate (decreases 
pressure) and the population of the thermal states (in¬ 
creases pressure). Moreover, this leads to a spurious lack 
of avoided crossing of the first- and second-sound veloci¬ 
ties, therefore necessitating the use of at least the Popov 
theory to describe the sound velocities accurately. 

The entropy per particle is presented in Fig. Since 
the system is described by a mixture of noninteracting 
phonon gases at low temperatures, the entropy per par¬ 
ticle shows a power-law behavior at low temperatures. 
Furthermore, the phonon velocity scales with the inter¬ 
action strength and therefore stronger-interacting gases 
have a higher entropy. 


III. SPIN DRAG IN A 
PARTIALLY-CONDENSED BOSE GAS 

In this section we investigate kinetic (incoherent) pro¬ 
cesses that contribute to density and velocity dynamics 
for the balanced binary Bose gas below the critical tem¬ 
perature. We start by generalizing the results of the ki¬ 
netic theory for a uniform single-species Bose gas (see, for 
instance. Ref. (621) to the two-species gas. The collision 
integrals that we consider below can be formally derived 
from the Heisenberg equations of motion for the atoms 
of the two different species. For a detailed discussion of 
that derivation for the case of a single species, we refer to 
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FIG. 1. (Color online) Pressure of the two-species balanced gas. The interaction parameter increases from left to right: 

= 0.01, 0.05, and 0.1. Intercomponent interactions are zero (7 = 0) for the solid line, moderate (7 = 1/2) for the dashed 
line, and strong (7 < 1 ) for the dotted line. 




FIG. 2. (Color online) The entropy per particle of the two-species balanced gas. The interaction parameter increases from left 
to right: r\}^^a = 0.01, 0.05, and 0.1. Intercomponent interactions are zero (7 = 0) for the solid line, moderate (7 = 1/2) for 
the dashed line, and strong (7 < 1) for the dotted line. A limited temperature range is depicted to emphasize the difference 
between the curves. 


Ref. [53]. The collision integrals discussed in this section 
conserve the number of f and | particles separately. 

Since the kinetic processes in question vanish for very 
low temperatures, and in order to avoid complications 
posed by the fact that the Bogoliubov transformation 
mixes the spin species, we work in the Hartree-Fock ap¬ 
proximation in this section. The Hartree-Fock approx¬ 
imation is therefore valid for the whole range of tem¬ 
peratures where the collision processes play a significant 
role. Moreover, here we consider the “local equilibrium” 
situation where the chemical potentials of the thermal 
atoms and the condensates are equal to fi = g+n + gn'. 
The thermal atoms thus feel the Hartree-Fock mean-field 
energy {g^ + g)n. These considerations lead to a Bose- 
Einstein distribution with a nonzero average momentum 
Pnct (Pnci) for the noncondensed f (i) atoms. 


■'* e'®(Pi“Pnc<^T/2m-l-,S(s+-l-s)"-/3M _ 1 ’ ^ ' 

where a labels the spin species and i labels the momen¬ 
tum variable (the omission of the i label denotes the 
momentum variable p). Furthermore, we allow for dif¬ 
ferent nonzero momenta of the two condensates: The f 
(4,) condensate has a momentum of p^^ (Pci)- 
In the single-species case only one collision process, 
namely C 12 , using the notation of Ref. |63|, is responsi¬ 
ble for the relaxation of the difference between the con¬ 
densate velocity and the velocity of the thermal atoms. 
However, five different processes (see Fig. can con¬ 


tribute to the dynamics of various densities and veloci¬ 
ties in the two-species system, leading to the following set 
of quantum kinetic equations for the distributions of the 
thermal particles, where we only consider the collisional 
contributions: 

aj^lcoii = ^2 + + ^^[2 + ^^2^ (62) 

Icon = (5^2^ + Cii + eg + Cil. (63) 

We now discuss the individual collision terms present 
in these equations. We use a shorthand notation for the 
momentum integrals f^ = J d^Pi in order to simplify the 
following formulas. First, we have a spin-drag term (Fig. 
1^ (a); cf. Refs. UH SHj), which only involves scattering 
between the thermal atoms and thus exists both below 
and above the critical temperature, 

C22 = [ M 2 [S^^Hp-P 4 )-d'^^\p-Pi)], ( 64 ) 

^1234 

^22 = [ ^22 (p - P3) - (P - P2)]: ( 65 ) 

-'1234 

where we have denoted the common part of the integrand 
by 

^22 = +P2-P3-P4) ( 66 ) 

y-Silpl + pj - pI - pl]/ 2 m) 

X + fl) - (1 + /g(l + g)/^/!] • 
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FIG. 3. Diagrams corresponding to the relevant collision integrals. Solid lines depict thermal atoms, while dashed lines represent 
condensed atoms. 


This term contributes to the relaxation of the difference 
of momenta p^ct ~ Pnci ■ Furthermore, we consider the 
following intraspecies collision terms [Figs. I^b) and||;c)], 
which explicitly depend on the condensate density and 
thus only exist below the critical temperature, 

^123 

- S’'^\p-P2) - S^^\p-P3)], (67) 

where 

^n(U) = + P-til) -P 2 -P 3 ) (68) 

x^([p? +Pcta) - P 2 - Pl]/‘^rn - gno^(i)) 

X - (1 + ■ 

These terms contribute to the dynamics of the conden¬ 
sate fraction and also describe the relaxation between the 
condensate velocity and the average thermal particle ve¬ 
locity of the same species. Moreover, we also have similar 
terms for the interspecies scattering, namely, 

c;««' = / [J«>(p - p,) - «BI(p - p,)], 

^123 

(69) 


C 


-n(tT) 

12 


J 123 


(70) 


where 


^ ^ -P 2 -P 3 ) (71) 

x<5([p? -b - P 2 - P3]/2w - 5no;(t)) 

X - (1 + . 


The latter terms also describe the dynamics of the con¬ 
densate fraction, as well as the relaxation of various ve¬ 
locities mediated by the condensate. 


In order to obtain the equations for the change of the 
density and the momentum of the thermal particles, we 


perform an integration of Eqs. (62) and (63) leading to 

(Pp 


= I 




(72) 

Pdtf^i)- (73) 


(27r/i)3 

d?p 

j {2'Kh)^‘ 

Adding the equations that change the thermal densities 
and the equations that change the condensate densities 
obtained from straightforward considerations of the col¬ 
lision integrals, we have 

dtn\ = rll + ftl, 

=rti + fli, 
a^not = -rll - f\l, 

dtnoi = - fll. 


(74) 

(75) 

(76) 

(77) 


where 


F 


F 


n(u) 

12 

mw 

12 


1 


(27rfi,)3 
1 / 


in(u) 

^12 I 


'123 


A 


n(n) 

12 j 


(78) 

(79) 


'123 


(27r/i)3 

and we see that the total densities n^(i) are conserved 
separately. In a similar manner, the equations for the 
change of momenta are 
1 

I /„ „ > /I 


(P4 -Pi)A 22 -PctFr 


/1234 


‘ + (80) 


(27r/i)3 

= -7t4i 3 / (P4 - Pl)^22 - Pc^F 

{ZTTn) J 1234 


ill 

12 


(27r/i)3 


+ (81) 


(27rli)3 

dtino^Pct) = Pct^ll + Pettit 
5t(notPct) = Pctri2 +Pctri2- 
Note that the total momentum is conserved; thus, 

5t(^tPnct + n[Pncl + mPet + noiPcl) = 0. 


(82) 

(83) 

(84) 


For the hydrodynamic theory in the subsequent sec¬ 
tion, we are interested in the linearization of the colli¬ 
sion integrals in terms of the velocity differences. It is 
















straightforward to show that all the integrals are at 
least quadratic in terms of the momenta and, therefore, 
in linear response the densities stay constant for both 
species and dtn'^ = dtUoa = 0. In addition, this result 
implies that the condensate momentum experiences no 
linear relaxation. This is consistent with the common 
physical intuition that the condensate motion should not 
decay in the lowest order, as the condensate motion cor¬ 
responds to the flow of a superfluid. 

25 

20 
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0 
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FIG. 4. (Color online) Condensate-assisted (solid line), ther¬ 
mal (dashed line), and total (dotted line) spin-drag rates for 
a system of sodium atoms at a density of 10 ^^ ni“®, 7 ~ 1 , 
= 0.03. 


Nevertheless, the momentum difference between the 
two species of thermal particles experiences relaxation. 


dtiPnct - Pnci) 


-( — + — ) (Pnct - Pnci) 

\T22 T 12 J 

-^(Pnct-Pnci)> (85) 

■7^_ _i ' ^ 


where we have also defined the total spin-drag rate 1 /Tgd 
in terms of the thermal and condensate-assisted spin-drag 
relaxation rates 1 /r 22 and 1 /ti 2 , respectively. While the 
thermal spin-drag relaxation time 


hp 

T22nka^^ 



1 1 f°° dqduiq^ 

(nA3)2 sinh^(a;/2) 

q‘^/ld7Z+^gno—uj/2+7TUJ^/q'^ _ 

^q'^/ 16 'K-\-^gno—uj/ 2 -\-TTLo‘^/q^ ^ 


2 


( 86 ) 


has been calculated before in Ref. [M] for the tem¬ 
peratures above the condensation temperature, the 
condensate-assisted spin-drag rate. 


hp 


64 no a 


Ti 2 nAa^^ 3(27r)3 n 

1 


a f 

- J dpi dp3 Pi pI 


X 1-h 


e(P?+Pi)/47r+2/33rao _ 1 

1 1 


ef 


PlP3 


gpf/4ir-|-/3gno _ X gPi/4ir-|-/3grio _ ^ V 27r 


(87) 



where 0 denotes the Heaviside step function, has not 
been investigated before. As we can expect, the ther¬ 
mal spin-drag rate l/r 22 dominates at high tempera¬ 
tures, while the condensate-assisted spin-drag rate l/ri 2 


is more important at low temperatures (see Fig.|^. Even 
though in our approximation the thermal spin-drag rate 
has a maximum at the critical temperature, a more care¬ 
ful calculation leads to its suppression due to critical 
fluctuations |65j in a very narrow temperature window 
around T^. Therefore, strictly speaking, 1 /ti 2 and l/r 22 
vanish at both T = 0 and T = T^.. 


IV. HYDRODYNAMIC THEORY FOR A 
BALANCED BINARY MIXTURE 

The goal of this section is to derive a set of hydro- 
dynamic equations for a balanced mixture of two compo¬ 
nents, where each of the components has a superfluid part 
and a normal part. Note that the following discussion is, 
in principle, not limited to weakly interacting Bose gases, 
since it is only based on conservation laws. However, we 
do not discuss various dissipative terms such as the ther¬ 
mal diffusivity or the (second) viscosity [T71 ISDl [Ml EZ] j 
with the exception of the spin-drag term. Furthermore, 
in the previous section we have concluded that the spin- 
drag term only relaxes the noncondensate spin currents, 
whereas in what follows we discuss the hydrodynamics in 
terms of superfluid and nonsuperfluid (normal) currents. 
Therefore, we posit that the nonsuperfluid spin currents 
relax in exactly the same manner as the noncondensate 
spin currents. This is an approximation, which is only 
valid in weakly interacting systems, where the conden¬ 
sate and the superfluid are very similar objects. 

In the equilibrium situation, the particle densities of 
the two components are identical for both the superfluid 
(n|^ = = rf^) and the normal fluid (n(l^ = = n"^), 

and therefore the total densities of each species are iden¬ 
tical, too (tt-x = n|^-|-n"^ = n^). Moreover, in equilibrium 

there are no particle currents = j’f = as 

all the velocities vanish: where 

the normal current of the f component is defined as 
and the other currents are defined simi¬ 
larly. The current of each of the components is the sum 
of the superfluid and normal currents: j-^ = + j’f ■ 

In the nonequilibrium situation, however, both the ve¬ 
locities and the various densities can be nonzero and dif¬ 
ferent from each other. In that case, it is useful to define 
variables pertaining to the combined motion of the whole 
gas and contrast them to the variables that describe the 
relative motion. Therefore, 

• nf.sf •nf,sf , -nfjSf /oo\ 

3 tot =3^ +3i , (88) 

are the total normal and superfluid currents of the com¬ 
bined motion, respectively. Adding these two currents up 
yields the total particle current + itot- Fur¬ 

thermore, we are now in a position to define the combined 
velocities such that 

(89) 

and 

ftit = «tot^tof 


(90) 
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Note also that thermodynamic functions (e.g., pres¬ 
sure and entropy) are dehned for the whole system only, 
and not for the individual components. Moreover, due to 
time-reversal invariance, the lowest-order velocity correc¬ 
tions to thermodynamic functions are quadratic |66j and 
therefore do not enter linear equations. For example, the 
Gibbs-Duhem relation reads 

-I- nid^i -I- sdT = dp (91) 

even in the case of nonzero velocities. Hence, we obtain 
the following linearized hydrodynamic equations for the 
combined (or in-phase) motion of the whole gas: 

dtntot + V • Jtot = 0, (92) 

dtiritots) + ntotsV ■ = 0, (93) 

+ Vp = 0, (94) 

+ Vptot = 0. (95) 

These simple equations express particle conservation, 
entropy conservation, Newton’s second law, and the 
Josephson relation, respectively. Note that the terms 
“in-phase” and “out-of-phase” in this paper refer to the 
motion of the two spin components, as opposed to the 
normal and the superfluid component, where the latter 
meaning is common in the literature concerning liquid 
helium. 


We now turn our attention to the relative motion of 
the up and down particles. To that end, we define the 


normal and superfluid spin densities 


An"*’**f = nf- nf 

(96) 

as well as spin currents 


Aj ’ =3^’ -3^' ■ 

(97) 


Using these newly defined quantities, the hydrodynamic 
spin equations are the following: 

-h V • = 0, (98) 


dtAn^^ + V • Af^ = 0, (99) 

dtAr^ + <ftVA/r/2m = -Aj^Vr^d, (100) 

dtAf^ + nli^VAp/2m = 0 . ( 101 ) 

The four equations above state that both normal and 
superfluid spin densities are conserved and the normal 
spin current is driven by a chemical potential difference 
and relaxes with the spin-drag rate, whereas the super¬ 
fluid spin current is also driven by a chemical potential 
difference but does not relax. 


By eliminating all the currents and velocities from the 
equations above and noticing that for the balanced case 
dAp 


we find that 


dT 


mdfS = 


dAfi^ 

dn 

(102) 

V'p, 

n"f ’ 

(103) 

(104) 


a?An"f -h —aAn"^ = 


'^sd 


(105) 


dAfj.^ 
m dAn 


V^(An''*-F An"f), 


9“ An"* = 


(106) 


dAfi 2f.nl 


V^(An'** + An"*). 

m oAn 

We now employ a traveling-wave ansatz for the total den¬ 
sity 


ntot = ntot,eq + Sntot exp(-*[a;t - k ■ x\), (107) 

where ntot,eq is the total density in the equilibrium, and 
Sritot is the amplitude of the wave. Proceeding similarly 
for the other quantities, we obtain the lowest-lying col¬ 
lective modes for the system. Note that in this balanced 
situation the in-phase collective modes (first and second 
sounds) are decoupled from the out-of-phase modes (spin 
modes). 


V. RESULTS 

In this section we present the results of our calcula¬ 
tions, most important of which are the sound velocities 
of the various modes. We discuss both the uniform and 
the trapped cases. From the total number of hydrody¬ 
namic equations, we expect four modes: two in-phase 
modes and two out-of-phase modes. The two in-phase 
modes are present in any system where entropy transport 
is distinct from density transport: They are the first- and 
the second-sound modes. It is worthwhile to note that 
the phenomenon of the second sound has been reported 
not only in superfluid helium, but in other systems as 
well, including solid helium [^, dielectric crystals |69) . 
and, more recently, a unitary Fermi gas M- Moreover, 
contrary to the case of superfluid helium, second sound 
in weakly interacting gases is not a pure entropy wave, 
as discussed below. On the other hand, the out-of-phase 
modes are similar to spin modes and therefore only occur 
in systems with several species of particles. 


A. Collective modes in a uniform gas 


In the case of no spin drag, all the collective modes of 
the system have a linear dispersion 

u! = Cik, (108) 

where ci are the sound velocities. At zero temperature, 
the first-sound velocity can be calculated from the ex¬ 
pression of pressure in Eq. (44), yielding 

o; = = <1±A5, (109) 

m m 

Moreover, the second-sound velocity can be calculated in 
a manner similar to that used in the single-component 
case (see Ref. for an explicit calculation), the only 
difference being that, instead of a single phonon gas, we 
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FIG. 5. (Color online) Square of the sound velocities of the two-species balanced Bose gas scaled to the zero-temperature 
second-sound velocity of a single species. Here = 0.03. The inter-component interactions are increasing from left to 

right: 7 = 0 (left), 7 = 1/2 (center), and 7 < 1 (right). The latter corresponds to a mixture of hyperfine states of the sodium 
atom with the density of lO^'^ m“^. The first-sound velocity is depicted by the thick solid line, the second-sound velocity is the 
dashed line. The spin-sound velocities are the dotted and the thin solid lines. 


have to consider a mixture of two noninteracting phonon 
gases, leading to 


c 


2 

2 


gn (1 -b 7)(1 - 7)^/^ + (1 - 7)(1 + 
3m (1 — 7)5/2-I-(1-I-7)®/^ 


When the two components are decoupled (7 = 0), we re¬ 
cover the single-species result c| = gn/3m from the for¬ 
mula above. Furthermore, in the strong-coupling limit 
(7 ^ 1)) sound velocity vanishes (c^ — 0), signal¬ 
ing the demixing transition. When it comes to the spin 
sounds, one of them always has a zero velocity, whereas 
the other one has the velocity 


g_^ 


(1 -7)gn 


( 111 ) 


as can be seen from the difference of chemical potentials 
in Eq. (22). Note that ci and Cg can also be obtained by 
expanding the energies of the Bogoliubov excitations in 
Eq. (551 to the lowest order in momentum. Therefore, the 
first sound and the spin sound can be thought of as both 
quasiparticles and collective modes at zero temperature. 


We present the velocities of the sound modes for differ¬ 
ent temperatures in Fig. In order to characterize the 
sound modes, we have calculated the amplitudes of the 
temperature and density perturbations from the eigen¬ 
vectors of the linearized system (see Fig. |^. We ob¬ 
serve that for the temperatures below the avoided cross¬ 
ing, first sound is mostly a temperature wave. However, 
above the avoided crossing temperature, first sound is 
predominantly a density wave. Second sound has compa¬ 
rable contributions from both relative temperature and 
density deviations for any temperature. Furthermore, 
second sound is a wave where density and temperature 
change in phase, while first sound describes an out-of- 
phase change (the temperature increases with decreasing 
density), as signified by a minus sign in Fig. When it 
comes to the spin modes, a zero-frequency mode exists, 
which corresponds to 

( 112 ) 

and does not affect the total density of either component. 
The other spin mode, however, affects the total spin den¬ 


sity, while the normal component contributes with the 
same relative weight as the superfluid component: 


An"f _ An"f 


(113) 


Upon accounting for spin drag, we find that the zero- 
frequency modes split into two: one zero-frequency mode 
and one purely imaginary mode. The purely imaginary 
mode (cf. Fig. at low momenta is an excitation with 
only normal density fluctuations, whereas at higher mo¬ 
menta it preserves the total density of every component 
as the zero-frequency mode: An“^ = — An®^. Further¬ 
more, the dispersion of the other spin mode develops a 
quadratic imaginary part, even though its real velocity 
Cg is almost unaffecte d by spin drag and the mode is still 
characterized by Eq. (1131. The frequency uj of the latter 


mode at long wavelengths can be written as 
w = Cgfc — iDk'^, 


(114) 


where k is the wavenumber and D is the diffusion coef¬ 
ficient. Hence, we call this mode the complex spin mode 
in Fig. 0 


Since all the interaction strengths are similar (7 ~ 1) 
in a sodium gas, when spin drag is absent, both spin 
modes have zero frequency. In particular, the spin sound 
has zero velocity at any temperature, Cg ~ 0, which is 
consistent with the zero-temperature result in Eq. ( 111[ ). 
When spin drag is present, the purely imaginary mode 
has a constant imaginary frequency ui = —ijTsd-, while 
the complex mode frequency remains zero, and no imag¬ 
inary part develops, so H ~ 0 in this case. In a re¬ 


alistic sodium gas with a density of 10^^ m”^ at half 
of the critical temperature, the total spin-drag rate is 
1/Tsd — 1 kHz, while in experiments of trapped sodium 
gas [48] above Tc typical rates are on the order of 0.1 kHz. 
Therefore, it should be experimentally possible to mea¬ 
sure the spin-drag rate below the critical temperature by 
measuring the spin-wave decay time. 
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FIG. 6 . (Color online) Ratio of density and temperature fluctuations present in first sound (solid line) and second sound 
(dashed line) for = 0.03. Intercomponent interactions are increasing from left to right: 7 = 0 (left), 7 = 1/2 (center), 

S'lid 7^1 (right)- The latter corresponds to a mixture of hyperfine states of the sodium atom with the density of 10^^ m“®. 



FIG. 7. (Color online) Spin modes in the — 0.03 case 

at T = Tc/2 for nonzero interspecies interactions (7 = 1/2, 
l/Tsd = 0.25kHz). The solid line depicts the imaginary mode, 
whereas the dashed and dotted lines depict the imaginary and 
the real parts of the complex spin mode. The horizontal axis 
corresponds to the physically relevant wavelengths down to 
about lOpim. 


velocities are presented in Fig. At zero temperature, 
the first-sound velocity in a trap is suppressed by a factor 
of \/2 as compared to the uniform case, since the aver¬ 
age density in the trap is half of the peak density in the 
Thomas-Fermi approximation, i.e., 

(no) = no/2, (115) 

as first pointed out by Zaremba m- Since at zero tem¬ 
perature the second-sound velocity is also proportional 
to the square root of the density, it is suppressed by the 
same factor of This is not explicit for sodium atoms 
as the second sound at zero temperature vanishes. How¬ 
ever, for a nonzero temperature we have not succeeded 
in finding a similar simple relation between the trapped 
and uniform gases. Finally, we also present the spin-drag 
rate dependence on temperature for a trapped system 
(Fig.§. The curves are qualitatively similar to the uni¬ 
form case, even though the rates are decreased by a factor 
of about 10. 


B. Gas in a trap 


We now turn our attention to the experimentally rele¬ 
vant trapped case. Here we consider a cylindrically sym¬ 
metric trap which is highly anisotropic. This trap, where 
one (axial) direction is very shallow, and the other two 
(radial) directions are much more strongly confined, puts 
the gas in the hydrodynamic regime in the axial direction 
even in the presence of weak interactions. Therefore, we 
perform the trap average in the radial direction and ana¬ 
lyze the propagation of excitations in the axial direction. 
To that end, we work in the semi-classical approximation 
(see Ref. [55] for more details). The condensate is treated 
in the Thomas-Fermi approximation and the excitations 
are treated in the Hartree-Fock approximation. There¬ 
fore, both components are treated in the local-density 
approximation. We consider the experimentally relevant 
situation m of 2-10® sodium atoms in a trap with an ax¬ 
ial trapping frequency of 27r x 2 Hz and a radial trapping 
frequency of 27r x 100 Hz, and extract the radial profile at 
the center of the trap from this calculation. We then per¬ 
form an average on all the thermodynamic quantities of 
the hydrodynamic Eqs. (103)-(106). The resulting sound 



FIG. 8 . (Golor online) Trap averaged sound velocities nor¬ 
malized to the peak density in the trap with no damping (see 
text for details). 


VI. CONCLUSION 

In summary, we have constructed a hydrodynamic the¬ 
ory of a balanced two-species Bose mixture. We have 
also calculated the microscopic thermodynamic parame- 
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FIG. 9. (Color online) Trap averaged condensate-assisted 
(solid line), thermal (dashed line), and total (dotted line) 
spin-drag rates normalized to the peak density in the trap 
(see text for details). 


ters entering that theory using the Popov approximation, 
obtaining the equation of state on the way. Moreover, we 


have accounted for the relaxation of the normal current 
by the spin-drag mechanism, considering the condensate- 
mediated spin-drag term for the first time. Adding these 
components together we were able to calculate the sound 
velocities and spin-drag rates for the experimentally ac¬ 
cessible system of sodium atoms in a trap. We hope that 
our analysis will stimulate experimental work on hydro- 
dynamic modes and spin drag in a partially condensed 
Bose mixture. 
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